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S_j ' 

We propose and analyze a scheme for the generation of multipartite entangled states in a system 

of inductively coupled Josephson flux qubits. The qubits have fixed eigenfrequencies during the 

whole process in order to minimize decoherence effects and their inductive coupling can be turned 

pO ■ on and off at will by tuning an external control flux. Within this framework, we will show that a 

W state in a system of three or more qubits can be generated by exploiting the sequential one by 

one coupling of the qubits with one of them playing the role of an entanglement mediator. 

PACS numbers: 03.67.Mn, 03.67.Lx, 85.25.Dq 
O 

o 

Vh ' I. INTRODUCTION 

Oh 1 
3 ■ 

Entanglement "the striking feature of quantum mechanics" (as claimed by E. Schrodinger in 1935 1 -) has been 

considered essential since the very beginning in order to investigate fundamental aspects of quantum theory. Quite 

recently, however, physicists have fully recognized that the generation of entangled states is an essential resource also 

in quantum communication and information processing. Entangled states have been generated in many experiments 

involving cavity QED and NMR systems, ion traps, and solid-state (superconducting) circuits, and their applications 

in the field of quantum computing have been demonstrated^ 3 *^*^ Among the previously mentioned physical systems, 

Josephson-j unction based devices presently provide one of the best qubit candidates for the realization of a quantum 

computer due to the fact that a wide variety of potential designs for qubits and their couplings are available and that 

qubits can be easily scaled to large arrays and integrated in electronic circuitsA& A series of successfully performed 

experiments for charge, flux, phase, and charge-flux qubits show indeed that they satisfy DiVincenzo's prescriptions^ 

^. , for quantum computing in terms of state preparation, state manipulation, and readout. Moreover, the nonlinearity 

characterizing Josephson junctions and the flexibility in circuit layout offer many possible options for coupling qubits 

i~~l , together and for calibrating and adjusting the qubit parameters over a wide range of values. 

In this field, remarkable achievements include the realization of complex single-qubit manipulation schemes^ the 
generation of entangled state s 11 ' 12 in systems of coupled flirsi 3 - and phas o 14 i 15 qubits, as well as the observation of 
quantum coherent oscillations and conditional gate operations using two coupled superconducting charge qubits>ii*i& 

(~ } The next major step toward building a Josephson-junction based quantum computer is therefore to experimentally 
realize simple quantum algorithms, such as the creation of an entangled state involving more than two coupled 
qubits, 3 . 6 . 17 ' 18 ! 19 

This goal may be achieved by selectively turning on and off the direct couplings between two qubits or their interac- 
tion with auxiliary systems (LC-oscillator modes^ inductances^ 1 - large-area current-biased Josephson junctions ; 22 ! 23 
or the quantized modes of a resonant microwave cavitjt 2 ^* 2 ^ 6 ) playing the role of a data bus. Typically, the coupling 
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energy may be controlled by tuning the qubit level spacings in and out of resonance. However, in order to avoid 
introducing extra decoherence with respect to that characterizing single-qubit operations, other promising scheme for 
realizing a tunable coupling of superconducting (spatially separated) qubits have emerged: for instance that wherein 
the interaction between two flux-based qubits is controlled by means of a superconducting transformer with variable 
flux transfer function^ or those wherein two qubits with an initial detuning can be made to (resonantly) interact by 
j^ ' applying a time-dependent (microwave) magnetic flux to the qubits^ 

Within these experimental frameworks, we propose a theoretical scheme by which it is possible to entangle more 
than two (spatially separated) flux qubits. It is based on the sequential inductive interaction of the qubits with one 
of them acting as an entanglement mediator. 29 ' 30 More in detail, we will see that the scheme operates in such a way 
that it generates an entangled W state after a finite number of steps, that no conditional measurement is required, 
and that the proposed architectures are scalable, at least in principle, to an arbitrary number of qubits. 

The paper is organized as follows. First, in Sec. ^ we briefly describe the main features and the Hamiltonians 
characterizing two kinds of Josephson devices, namely the double rf-SQUID 3 ^ 3 ^ 3 " and the persistent (three-junction) 
SQUID, 34 by which it is possible to build superconducting flux qubits. Then, we discuss the most common experi- 
mental procedures by which it is possible to initialize and to measure their quantum state. In Sec. IIIII a scheme of 



successive interaction is introduced in an (iV + l)-qubit system, qubit M+qubit l+< ■ ■ +qubit N, wherein entanglement 
mediator M is coupled one by one with qubits 1, 2, . . . , N. We discuss moreover the possibility of practically realizing 
this scheme by exploiting some of the different physical coupling elements currently available and how the coupling 
energy depends on the particular way in which the interaction between each qubit and the mediator is implemented. 
In Sec. IIVI we analyze the dynamics of the system showing that, by preparing the multi-qubit system in a pure 
factorized state and by adjusting the coupling energies and/or the time of interaction between each of them and the 
mediator, an entangled W state can be generated. Finally, conclusions and discussions are given in Sec. 

II. SUPERCONDUCTING FLUX QUBITS: MODELS AND HAMILTONIANS 

In this section, we briefly describe the main features and the Hamiltonians of two devices, the tunable rf- 
SQUID2ii22i2£ and the three-junction SQUID, 34 by which it is possible to implement a two-state Josephson system, 
focusing our attention also on the experimental procedures to be considered in order to prepare their initial quantum 
state. 

A. The double rf-SQUID qubit 

We begin by considering a double rf-SQUID system; 3 - 1 * 3 ^ 3 - that is a superconducting ring of self-inductance L 
interrupted by a dc-SQUID, a smallest loop containing two identical Josephson junctions, each with critical current 
io and capacitance c. This device [schematically illustrated in Fig. ^a)] is biased by two magnetic fluxes 4>x and 4> c 
threading the greatest ring and the dc-SQUID, respectively. 

The dc-SQUID if small enough (i.e. with inductance I <C iq4>q/2tt 1 where <p — h/2e is the flux quantum) behaves 
like a single junction with a flux-dependent critical current I c — 2io\ cos(7n^ c /<^o)| and capacitance C = 2c. This 
means that a double SQUID simulates a standard rf-SQUID with tunable critical current I c = I c (<fi c ). 

Therefore, by taking into account both the charging energy of the "effective dc-SQUID junction" (T = q 2 /2C) and 
the washboard potential, the Hamiltonian of the system is written down as 

H= 2C + ^2L ^T COS Uo-J' (1) 

where the charge on the junction capacitance, q, and the flux through the SQUID loop, </>, are canonically conjugate 
operators satisfying the commutation relation [<fi, q] = ih. 

It is well known 33,35 that, by setting /3l = 2irLI c /(f)o > 1 and <fi x « </>o/2, the circuit behaves as an artificial quantum 
two- level atom whose reduced Hamiltonian in the basis of the flux eigenstates \L) and \R) (which are localized in the 
two minima of the washboard potential and correspond to two different orientations of the current circulating in the 
large loop) assumes the form 

Hrf = ~-zA KF ((j) c )a x - -eRF(0x)cr z . (2) 



Here, he-R.v{<t>x) = 2I cy /6(/3i J — 1)(4> X — 4>o/2) is the energy difference between the two minima of the washboard 
potential, 

(in the limit < (3l — 1 C 1) the tunnelling frequency between the left and the right wells that can be tuned by 
changing the junction critical current I c ((j) c ), and a x and a z the Pauli spin operators. 

B. The persistent-current (3 J J) qubit 

To minimize the susceptibility to external noise of a large-inductance rf-SQUID, Mooij et al£^ proposed to use a 
persistent-current qubit [schematically shown in Fig.^b)], namely a smaller superconducting loop containing three 
Josephson junctions, two of equal size (i.e. with Eji = Ej 2 = Ej) and the third smaller by a factor a (i.e. with 
-^J,3 = a Ej, a < 1). By applying an external magnetic flux 4> x close to </>o/2 and choosing a ~ 0.8, it has been proved 
that, in the low-inductance limit (in which the total flux coincides with the external flux and fluxoid quantization 
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FIG. 1: Sketch (a) of the superconducting double rf-SQUID qubit and (b) of the three-j unction SQUID by which it is possible 
to realize a persistent-current qubit. 

around the loop imposes the constraint ip\ — </?2 + <^3 + 27r/ = on the phase drops across the three junctions, 
/ = fix/fio being the reduced magnet flux), the Josephson energy 



U(tpi,(p 2 ) = -Ejcosipi - Ej cosip 2 - aEj cos(2Trf + ipi — tp 2 ) 



(4) 



forms a double well which permits two stable configurations of minimum energy corresponding to two persistent 
currents of opposite sign in the loop. This means that, also in this case, we may engineer a two-state quantum system 
whose effective Hamiltonian, in the basis of the two states that carry an average persistent current ±J p ss ±2naEj/(j)o 
(named \L) and \R) also in this case), reads 



TT h A ft I ± \ 



(5) 



Here, the tunnelling matrix element between the two basis states, ftA 3 jj/2, depends on the system parameters, and 
£3Jj(0x) = 2Ip(<f>x - <f>o/2)/ft- 



Initialization and readout 



Using as a new basis that is spanned by the energy eigenstates |0) and |1) (which are the symmetric and the 
antisymmetric linear superpositions of \L) and \R), if <j> x is exactly equal to </>o/2), the Hamiltonians of both the 
rf-SQUID qubit and the 3JJ qubit take the diagonal spin 1/2-like form 
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-ujo z 



(6) 



Here, hw = E\ — E — hy/e 2 + A 2 indicates the energy separation between their corresponding eigenstates (the 
analytic form of e and A, as previously discussed, depends on the specific design of the qubit) and <j z = |1) (1| — |0) (0| 
is a Pauli operator. 

Before discussing the modus operandi of our scheme, it is worth noting that both the rf-SQUID qubit and the 3JJ 
qubit are easily addressed and measured. Usually, they are initialized to the ground state simply by allowing them 
to relax, so that the thermal population of their excited states can be neglected. The coherent control of the qubit 
state is instead achieved via NMR-like manipulation techniques, i.e. by applying resonant microwave pulses which, 
by opportunely choosing the interaction time and the microwave amplitude, can induce a transition between the two 
qubit energy levels^ 

In addition, a flux state can be prepared producing a collapse of the wave function through a flux measurement or, 
as recently pointed out by Chiarello^ in the case of a double rf-SQUID qubit with an opportunely chosen sequence 



of variations of the washboard potential. A double SQUID can be indeed prepared in a particular flux state by 
strongly unbalancing the washboard potential in order to have just one minimum, then waiting a time sufficient for 
the relaxation to this minimum and finally sufficiently raising the barrier in order to freeze the qubit in this state. 
Finally, coherent rotation between the two flux states can be realized by lowering the barrier in order to induce fast 
free oscillations, waiting for fractions of the oscillation period to realize the desired rotation and opportunely raising 
the barrier in such a way to freeze the system in the desired target state. 

Also for qubit readout, several detectors have been experimentally investigated. Most of them include a dc-SQUID 
magnetometer inductively coupled to the qubit to be measured by which it is possible to detect the magnetization 
signal produced by the persistent currents flowing through it, exploiting the fact that the dc-SQUID critical current 
is a periodic function of the magnetic flux threading its loop^ In addition, it is worth emphasizing that, besides 
these proposals, physicists have been focusing their efforts on the realization of nondcmolition-measurement schemes 
(necessary for applications where low back-action is required) like that based on the dispersive measurement of 
the qubit state by coupling the qubit nonresonantly to a transmission-line resonator and probing the transmission 
spectrumi2£ 

III. THE SCHEME FOR ENTANGLEMENT GENERATION: SEQUENTIAL INTERACTION OF N 
FLUX QUBITS WITH AN ENTANGLEMENT MEDIATOR 

In this section, we propose a scheme for the generation of maximally entangled states in a multi-qubit system, 
M + 1 + 2 + • • • + N, where qubit M, playing the role of an entanglement mediator, is assumed to interact one by one 
with 1,2, . . . , and N. Among the different forms of coupling theoretically proposed and experimentally realized, we 
consider those by which it is possible to realize an inductive interaction2L22i2i between each qubit and the mediator, 
so that the free Hamiltonian of the whole system and the interaction Hamiltonians can be cast (in the basis of the 
energy eigenstates) in the following form: 

Ho= £ \^f (7) 

i=M,l,2,...,iV 

and 

H' Ml = 9l oi M) ^\ H> M2 =g 2 aMa<£\ ..., H' MN = g N a^a^. (8) 

We assume moreover that these qubits are properly initialized, by exploiting one of the previously mentioned ex- 
perimental recipes, so that at t = the whole system is described by the pure factorized state |1m0i02 . . -Oat) = 

|1)m<8> |0)i<g> |0) 2 <8> • • • <8> |0)jv 

With this setup, if qubits 1,2,..., and N are spatially separated (in order to strongly reduce their direct persistent 
coupling) and their interaction with mediator M can be turned on and off at will (by adjusting the coupling energies 
<7i, (72, • ■ ■ , 9n), we realize a step by step scheme which is sketched as follows: 

• Mediator M prepared in the state |1)m interacts inductively (by setting gi ^ 0) with qubit 1 during 
an opportunely chosen interval of time < t < T\, while qubits 2, 3, ..., and N evolve freely (with 
.92 = 53 = • ■ • = 9N = 0). 

• At t = n, we turn off the interaction between M and 1 and we adjust g^ in order to couple the mediator and 
qubit 2 (by choosing g^ ^ 0) for T\ < t < t\ + r%. 

• In a similar manner, we put qubits 3, 4, . . . , and N in interaction with M one by one. 

• Finally, at t = T\ +T2 + • • ■ + tjv, we switch off the interaction between qubit N and mediator M, and the desired 
entangled state of the (N + l)-partite qubit M+qubit 1+- • ■ +qubit N system is generated, provided that the 
interaction times, Tj (with j = 1,2, . . . , N), and/or the coupling constants, gj, are accurately selected. 

At this stage, it is important to consider a realistic experimental setup by which it is possible to implement our 
scheme. We begin by considering an inductive mediator-qubit coupling realized by means of a superconducting 
transformer with variable flux transfer function as in the paper of Cosmelli et aim* They show indeed that, by using 
a superconducting flux transformer modified with the insertion of a small dc-SQUID, it is possible to control the flux 
transfer function and therefore the inductive coupling constant, by modulating via an externally applied magnetic 
flux cj) cx the critical current of the dc-SQUID (see Fig. [2J . More in detail, they prove that the transformer can operate 
between two states with very different behavior: the "off" state where the transfer ratio 1Z = 4> ut/4>\n is minimum 
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FIG. 2: Sketch of a controllable flux coupling (CFC) circuit. A signal flux is applied on the left side of the CFC and its response 
is read out/coupled to the SQUID on the right. The control of the transmitted flux <j> ou t is achieved by modulating the flux 
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FIG. 3: Sketch of the scheme for the controlled coupling of more flux qubits proposed by Cosmelli et alw— 



(~ 0-0.1) and the "on" state with a transfer function ratio which may be also larger than 1. Under such conditions, 
it is possible to conceive an experimental scheme, like that depicted in Fig. [31 where the coupling between mediator 
M and the jth qubit may be effectively turned on by adjusting the control fluxes of the relative "switches" in such a 
way that one obtains gj ^- with all the other "switches" kept in the off state. 

Similarly, Plourde et al. propose to adjust the coupling strength characterizing the interaction between two 3JJ flux 
qubits by changing the critical current of a dc-SQUID (which is coupled to each of these two qubits) , finding that 
their coupling constant can be changed continuously from positive to negative values (and enables cancellation of the 
direct mutual inductive coupling if the two qubits are not spatially separated)^ 

Adopting one of these experimental coupling setups, it is therefore possible to realize the step by step scheme in 
the system of N + 1 qubits whose Hamiltonian during each of the aforementioned steps assumes the following form: 



Hmj — Ho 



H Mj - H M 



+ H 1 + H 2 + 



H 



N 



H Mj 



(j= 1,2,. ..,7V). 



(9) 



We observe that the structure of the multipartite Hamiltonian during each step allows us to simplify the study of its 
dynamics by confining ourselves to the analysis of the dynamics of the bipartite mediator- S QUID- 1 system during 
the time interval < t < t\ , to that of the mediator-SQUID-2 system during the second period n < t < t\ + t 2 , and 
so on (of course, provided that the free evolution of the other qubits is carefully taken into account). 

Moreover, by assuming that all the qubits and entanglement mediator M have a common energy gap ui = w,-, 
V i = {M, 1,2,..., TV}, due to the preponderance of rotating- wave terms of the interaction Hamiltonians with respect 



to the counter-rotating ones, it is not difficult to convince oneself that, during each step, the system dynamics is 
dominated by the bipartite Hamiltonian 

H^ A = H M + H t + H^ A ' = ^ + \»i<& + 9j (ai M) a^ + <T™<rf), ,m = v, = u,, (10) 

describing the rotating-wave coupling between the mediator and the jth qubit. In the following section, we will 
demonstrate the validity of this assumption. It is interesting, however, to note that Hamiltonian i|10|) can be exactly 
implemented if the detuning A = u>m ~ ^j between the two flux (3JJ) qubits to be coupled is chosen sufficiently large 
so that initially each of them can be treated independently. As recently shown by Liu et al.f& in fact, this gap can 
be nullified by applying to one of the two SQUID loops a time dependent magnetic flux 4> x (t) = Apcosuj c 't (with 
£ = M or j) satisfying the condition A ± lu c ' =0 when A ^ 0. This means that, by considering the reduced bias 
flux of each qubit close to (f>o/2 and the aforementioned frequency-matching condition, it is possible to implement 

the interaction Hamiltonian H 1 ^ = gj(<r + a_ + a_ <t_j? ), where gj can be controlled for instance by tuning the 
amplitude A( of the applied time-dependent magnetic flux. 

IV. W-STATE GENERATION 

Now, we analyze the dynamics of the system and demonstrate that the scheme introduced in the previous section 
actually generates an entangled W state of a tripartite system. We show further that it can be extended to the case 
of a larger number of qubits. 

We begin our analysis by looking at the eigensolutions of Hamiltonian (|1(J|) . which is represented in the basis 
{|01)mj, |10)mj, |00)a/j, |11)mj} by the following simple diagonal block form: 

/ o 9j 

f/EWA _ [ 9j , 

Hm ? - -fej I • ! " ' 

\ faj 

These blocks describe three dynamically separate subspaces: the first with characteristic frequency gj/fr characterizing 
the appearance of the entanglement between the degenerate states \01)mj and |10)j\/j, and the other ones describing 
the fact that the two states \00}mj and |11)mj evolve freely. We easily find that the eigenstates of Hfff A are 

\u lj ) = -j=[\10) Mj -\01) Mj ], \u 2j ) = ^=[\10)m j + \01)m j ], K-) = |00) Mi , \uaj) = \U) Mj , (12) 

with eigenvalues 

^ij/2j = Tffj, ^3j/ij = T^J- (13) 

A. Generation of entangled W states of the tripartite "M+l+2" system 

By exploiting the knowledge of the eigensolutions of Hamiltonian H^j A , it is possible to follow step by step the 
dynamics of the three-qubit system characterized by the one by one interaction of mediator M with qubits 1 and 2. 
We choose as an initial condition the state 

|1m0i0 2 ). (14) 

During the first step, we switch on the interaction between M and 1 and we let them interact for a time n, while 2 
evolves freely. This process generates the state 

e -«**m/ft| lAf0l02 ) = e -*(»o+ffg? A ')n/ft| lilf0l0a ) = e ^/ 2 [cos^|l M 0i0 2 ) - ism^|0 M li0 2 )], (15) 

where 9\ = g\Ti/fi. Next, by turning on the interaction between M and 2 at t = T\ and by allowing qubit 1 evolves 
freely during the second step (i.e. during the interval of time t\ < t < t\ + r 2 ), we obtain: 

\ V2 ) = e - 4iW2 / ri e-^" ri / R |l M 0i0 2 } 



_ iuj(T 1 +T 2 )/2 



[cos0icos6> 2 |1mOiO 2 ) - isin6»i|0 M li0 2 ) - icos6»i sin6» 2 |0 M 0il 2 )]. (16) 



Equation 116(1 clearly shows that, by adjusting 8j = gjTj/h (j — 1, 2) so that 

|sin#i| = ^P and | cos ^2 1 = —7= , (17) 

v 3 v 2 

a tripartite W state is generated. If we choose 9 2 — 7r/4 and d\ such that sin#! = l/\/3 and cos6>i = \/2/3, for 
instance, we get 

\W 2 ) = ^=e iu ^ +T ^ 2 [\l M 0i0 2 ) - i|0 M li0 2 ) - i|0 M 0il 2 }]. (18) 

It is worth noting that we obtain this state by adjusting the coupling energies gj during the aforementioned steps 
and/or by tuning the interaction times Tj. 

B. Generation of entangled W states of the multipartite "M+1+2+- • • -\-N" system 

Within this framework, it is possible to look at the possibility of applying the same techniques in order to generate 
an entangled state of a multipartite "M+1+2+- • • +A" system with N > 2. As described in Sec. Ill, we consider an 
entanglement mediator M in interaction one by one with TV qubits. If the system is prepared at t = in the factorized 
state 

|1 m 0i0 2 ...0jv), (19) 

after a straightforward approach it is easy to show that at the end of the iVth step (namely at t = ijv = X}i=i T j) ^ 
can be described in terms of the state 

\<Pn) = e^" 1 ^*"/ 2 [cos 6»i cos 6i 2 • • •cos6» A r|l A/ 0i0 2 ...0 N ) 

— «sin#i|0Mli02 • • -On) — • • • — icos6>i cos #2 • ■ -cos^-i sin6>fc|OjvfOifJ2 . . . 1& . . . Oat) 

i cos 9\ cos 6*2 • • • cos 9n-i sin 6>jv |0m0i02 ■ • • 1^)] • (20) 

Assuming that it is possible to control the interaction time tj of each qubit with mediator M and/or their coupling 
constants gj so that 



1 N-j+1 

sindi = and cos#, = \ , (21) 

J y/N-j + 2 3 p-J+2 ' 

we finally find that the generalized W entangled state 

i(N-l)ut N /2 

\W N ) = — [|1m0i0 2 ... Oat) - i|0 M li0 2 ...0 N }- i\0 M 0ih ■■■On) i|0 M 0i0 2 . . . 1^)] (22) 

V iv + 1 

of the (N + l)-partite system is created. 

C. Estimation of the effect of counter-rotating terms 

At this stage, we wish to test the validity of the rotating-wave approximation (RWA) performed at the end of Sec. 
IIIII To this end, we analyze the fidelity of the system state | V'jv) calculated without performing the RWA on the 
Hamiltonian model with respect to the target state (12211 . i.e. Tm = \(Wm\'4 ) n)\- 

Let us first look at the fidelity for the tripartite W state. After a straightforward calculation, we find that, by 
following the aforementioned two-step procedure, the a x oi couplings generate the state 



|^2) = i= e -^+^)/2 J |1 m0i 2 ) - z|0m0i1 2 ) - ie~ 1 ^ 



cosx2 + z sinx2 |0 M li0 2 ) 

v y/gi + t?LU 2 ) 

-t , 9 ^ ^r sinx2 ImIiIj 



(23) 



8 

with the same 6\ as that for i|18|) . where \2 — \Zg 2 /^ 2 + <jj2 ' T 2- Its fidelity Ti — 0.999957, calculated with gijh ~ 
0.5 GHz, wkIO GHz (in agreement with the currently available experimental values), confirms that during each step 
the system dynamics is dominated by the bipartite Hamiltonian ((THf describing the rotating-wave coupling between 
the mediator and the jth qubit. 

For the (N + l)-partite system, the generated state reads 



l<M 



i{N-l)uit N /2 



y/N + T 



|ljjf0 1 2 ...(M-i|0 Jtf 0i02...1 



N) 



N-l 



N 



-<e n «- 

fc=l \j=k+l 



Tj cos xj + i 



huj 



g 2 + h 2 Lu 2 



sm.Xj } |0M0i0 2 ...lfc...0jv) 



(states orthogonal to |Ww)) 



(24) 



with the tuning i|21|) . where \j — \/ y 2 /^ 2 + < jj2 ' T j-> an d the fidelity is given by 



F> 



N 



1 

N+l 



N-l N 

2+e n 



e lUT] I cos Xj + 1 



^g 2 + tfu 2 



'Xj 



which results in 

T z = 0.999873, T A = 0.999601, 



J"m = 0.997561, 



T20 = 0.993959, etc. 



(25) 



(26) 



when gi/h = ■ ■ ■ = gjq j% ~ 0.5 GHz and lo w 10 GHz. The physical meaning of such a list of values is that we might 
implement our scheme using up to 20 qubits maintaining a very good level of fidelity of the state given by Eq. I|24l) 
with respect to the target state expressed by Eq. I|22|) . 

V. CONCLUSIONS 

In this paper, we have discussed a scheme for the generation of entangled W states among three Josephson (eventu- 
ally spatially separated) flux qubits as well as its generalization to the case of TV + 1 qubits. The success of the scheme 
relies on the possibility of realizing controllable couplings between the qubits and entanglement mediator M and on 
the possibility of preparing their initial quantum state. It should be stressed that no conditional measurement are 
required but we have to tune the coupling energy and/or the interaction time between each qubit and the mediator. 

Final considerations are devoted to the analysis of some aspects of decoherence problem in our system. In particular, 
we wish to verify if the operation duration (e.g. the time necessary to perform the desired quantum process) is small 
enough with respect to the decoherence time. We find, on the one hand, that the eigenfrequency of a Josephson qubit 
lj is of the order of 10 GHz and that correspondingly Qj/U ~ 0.5 GHz. Under such conditions, the length of each step 
(during which only a fraction of a Rabi oscillation takes place) is at most of the order of h/gj as 2 ns and consequently 
the whole process in the case, for instance, of a tripartite system lasts approximatively 4 ns. On the other hand 
the relaxation and decoherence times, T\ and T2, of a superconducting flux qubit have been coarsely estimated and 
are in the range 1-10 fj,s&^ Therefore this means that a number of operations can be performed before decoherence 
occurs and that, in principle, our quantum information processing can be exploited for the experimental generation 
of entangled W states among more than three qubits. 
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